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Abstract 

The theory of stochastic differential equations is exploited to derive the Hopf's 
identities for a passive scalar advected by a Gaussian drift field delta-correlated in 
time. The result holds true both for compressible and incompressible velocity field. 
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1 Introduction 



In the last years much attention (see for example [|T], |2|, ^ ||, [|] and references therein) has 
been devoted to the study of passive scalar advected by an external field v = v(x, t) and 
described by the equation 

d t T+(v-V)T-vAT = f (1) 

where v denotes the molecular diffusivity of the scalar T = T(x, t) and / = /(x, t) is an 
external source driving the system. The solution of ([!]) in a time interval [0,t] is easily 
computed once it is known the kernel of the semigroup operator defined by the backward 
Fokker-Planck equation 

d t P + V ■ (vP) + uAP = 

limP(y, S |x,t) = 5(x-y) (2) 

s\t 

describing the conditional probability density for a particle to be in position y at time 
s < t when it arrives in x at time t. 

The connection between the two equations is provided by the trajectories describing 
the realizations of a stochastic differential equation backward in time (BSDE), which 
represents the Lagrangian picture of the passive scalar equation (section Q). This is the 
generalisation of the fact that in order to solve equation ([J) for zero viscosity one has to 
invert the solutions of the ODE generated by the flow v(x, t) . 

The general solution of (|2]) can be written in the form of a path integral. This obser- 
vation is very useful in the case of a random velocity (drift) field Gaussian in space and 
delta-correlated in time (section |3|). There it is shown that for each N-points function of 
the theory i.e. for each average over the drift of the products of N fields at equal times, 
the Gaussian average defines a system of BSDEs of N variables. The equations satisfied 
by the N-points functions are the backward Kolmogorov equations corresponding to each 
of the BSDEs systems (section |]). 

2 The Lagrangian picture 

Let us consider the stochastic differential equation 

dx s = — v(x s , t — s)ds + \f2~udw s s < t 
x = x (3) 

In general the solution of (||) exists and is unique if the drift field is Lipschitz and does 
not grow at infinity faster than |x| 2 ([§], 0). It describes the motion of a particle localised 
at time t in x and diffusing for increasing s i.e. backward in the global time s' = t — s 
From the rules of Ito calculus the total differential of a scalar functional of (g|) is given 

by: 

d s T(x s ,t- s) = 

= -[d t - s + v(x s ,t- s) ■ V - vA]T(x s ,t- s)ds + V2vdw s - VT(x s ,t-s) (4) 
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The average of a Ito stochastic integral is zero therefore if T satisfies ([I]) with a given 
initial condition T (x) at time s = its expression at time t is provided by Q: 

T(x,t)= / d D yT (y)P(y,0\x,t)+ f ds f d D yf{y,s)P(y,s\x,t) (5) 

JfL JO JR 

where P is the transition probability density specified by (§): 

P(y,s|x,t) =< 5(y-x t _ s ) > (6) 

If one uses again Ito calculus 

d s P(y, s | x, t) =< - d t - s 8(y - x t _ s ) >= 
= < [v(x*_ s , s) • V xt _ s - ^A Xt _J5(y - x*_ s )] > (7) 

by exploiting the properties of the delta function one arrives to (§) as mentioned in the 
introduction. Therefore the solution to the problem of the passive scalar is equivalent to 
the investigation of the problem of backward diffusion. 

The transition probability @ can be expressed in the form of a path integral: 

P(y, s | x, t) = (S) f t=X Vx u exp[- f ( l|x " + v <*» u ^ - - ly • v(x u , it) ) du] (8) 

The stochastic integral (§) is evaluated according to the Stratonovich's mid-point pre- 
scription as stressed by the symbol (S). In the case of additive noise the Stratonovich 
prescription has the advantage to preserve under path integral sign the rules of ordinary 
calculus. 

As it will become clear in the following, it is more useful for the scopes of the present 
paper to resort to the Ito pre-point prescription which allows to reabsorbe in the measure 
the divergence term in the short time Lagrangian of ([8]): 

P(y, a | x, t) = (J) / Vx u exp[- du\ (9) 

ix t _ s =y Jt-s 4 Z/ 

3 The Gaussian model 

Let us now assume that the velocity field in ([[]) is a Gaussian random field with zero 
average and correlation: 

< u a (x, fMy, s) >= 5(t - s)D Q ^(|x - y|) (10) 

Similar hypothesis are done for the external forcing: 

</(x,i)/(y, S )>=«5(t- S )C(x-y) (11) 

The assumption of time decorrelation is crucial in order to preserve under Gaussian aver- 
age the locality in time of the short time action in On the other hand no assumption 
is done on the incompressibility of the velocity field in contrast with . 
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In order to compute the general N-points function one needs to evaluate the path 
integral: 



<n^ 1 P(y,,0|x,,t,) > v = 
= f Dv(x, u) expf-S^II^ r ( * i)_X ' V^VX^u) exp[-S(x h X t )} (12) 



where 

rtMax f n d — y\)5(t — s) 

^( v ) = / duds Id xd yv a Cx,u) Wy, s) (13) 

Jo Jr 2 



and, after the Hubbart-Stratonovich transform: 

ft Max ... 2 



S(x,,A,) = / du0(t, -u)[i/|A f («)r + «Ai(u) -±i(u) + iXi(u) ■ v(x i («),«)] (14) 





In (|I4"D tjvjfoa: coincides with the largest ti and 9(ti — u) is the usual step function. The 
drift field appears linearly in (0): the Gaussian average can be easily performed. The 
averaged action is given by 

= T, l=1 du{X lta {u) + ^ Xi >a {u)xi >a {u)\ + 

Jo 2 

rmin( tl ,t k ) D a/3 (\ Xl (u) - X fc (n) |) 

+ Z l=lk>l / dwA iia (w) A fc)/3 (w) (15) 

where the Einstein convention and Greek letters have been used for vector indices (a, (3 = 
1, ...D). A further simplification is introduced by setting t\ = tW I: 

s av ({xi}l 1 ,{ML 1 ) = 

1 du {*M G ^ Mu l~^ u)) *Au) + * \Au)xiAu)\ (i6) 

with the Einstein convention now extended to all indices and 

G ajP (xi - x fc ) = G a>/3 (0) = [2 z/5 Q/3 + D aP (0)] I = k 

G Qi/3 (x; - Xjfe) = A^flxj - Xjfel) / 7^ & (17) 

From (|16l), (|T7|) is immediately clear that the one point function is Gaussian with constant 
diffusion matrix given by G ai p(0). On the other hand non-trivial behaviour is expected 
for N > 2. By integrating over the ghost trajectories {A;}^ one obtains 

Q IS \* \ I' A ■ I ^ G «j)(Mu) ~ Xk(u)) ■ ( s / 1a x 

S av ({xi\ l=1 ) = / duxi ja {u) — - xi,p{ u ) ( 18 ) 

Jo 2 

The averaged action (18) is exactly the one that we would have found starting from the 
system of Ito stochastic differential equations (SDE): 

dx lja (t) = Glpfait) - x k {t))dw ktP (t) l,k=l,N (19) 
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where dw(t) is the Wiener differential. For N=l it is again evident that the diffusion is 
normal with diffusion matrix G a ^(0). The effect of the Gaussian average is therefore to 
associate to each of the N points function an effective diffusion. 

Two remarks are needed. First, the system ( |19"D is autonomous. This implies that the 
relation between the transition probability density for the forward process and that one 
of the backward process is particularly simple: 

= ^ N) ({xi}lu{yi}li,t-s) (20) 

for any t > s since that 

^ {N) (Mi=i,{yi} = n^(x,- yi ) (21) 

Second, had we used the Stratonovich prescription we would have got to the averaged 
action: 

s av ({xi}D = 

= / du[xi ja (u) x k) p(u) - -di a d k pD a p(\xi(u) -x k (u)\) (22) 

Jo 2 4 

Whilst its probabilistic meaning is a posteriori clear, fl2"2| ) has the disadvantage of a less 
intuitive form. 

4 The Hopf 's identities 

Let us define the average of a given function F (x, ...,x) over the realizations of flTD| ) 
interpreted as BSDEs: 

< F( Xl (t), x N (t)) >= Uf =1 f°° d D yi F {x, x)pf >( yi , ...,y N ,s\ x 1; Xjv , t) (23) 

J— oo 

Then 

8 t < F(x x (t), ..., Xjv (t)) >= ^G a ^ - x k )d ha d k , p < F(x t \ ...,xf ) > (24) 

where |G Qi/ g(x/ — x. k )di >a d k) /3 is equal to the operator — M.j^ of f|, [|, |j, || Equation (|24]) 
is the backward Kolmogorov equation for the system (|1~9| ) interpreted as BSDEs so it is 
perfectly consistent with initial conditions in time. 

The Hopf's identities are a straightforward consequence of this observation. Namely 
starting from (H), let us consider the case of the two point functions. 

8 t «T(x 1 ,t)T(x 2 ,t) > />v = 
= fds r d D y 1 d D y 2 [T (y 1 )T (y 2 )5(s) + C(y 1 -y 2 )]d t V^\y 1 ,y 2 ,s\^ 1 ,^ 2 ,t) + 

JO J-oo 

+ d D y 1 d D y 2 C(y 1 -y 2 )vi > 2 \y 1 ,y 2 ,t\x u x 2 ,t) (25) 

J — oo 
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where < ... >/ denotes the average over the external force. By means of (|24|), we recognise 
that (H) can be rewritten as: 

d t «T(x 1 ,t)T(x 2 ,t) > f > v =-M 2 «T(x 1 ,t)T(x 2 ,t) > />v +C(x 1 -x 2 ) (26) 

In the limit t | oo (stationary case) we obtain the Hopf's identity for the two point 
function. 

Let us turn to the four point function. 

<9 t «nf =1 T(x / ,t) >/> v >= 

rt roo 

= dsdu Ut =l d D yi [nUT ( yi )S(s)8(s -u) + T ( yi )T (y 2 )5(s)C(y 3 - y 4 ) + 

JO J-oc 

+ C(yi - y2)C(y 3 - y 4 ) ] d t V i b 4) ({y 1 ,y 2 , s}{y 3 , y 4 , u] | x x , x 4 , t) + perm. + 

rt roo 

+ C(xi-x 2 ) / ds d D y 3 d D y4T (y 3 )T (y 4 )5(s) + 

JO J-oc 
(2) 

+ C(y 3 - y^lVb (Y3,y4,s | x 3 ,x 4 ,t) + 

+ perm. (27) 

here perm, means permutation among the / — indixes. 

The expression V b 4 \{yi,y 2 , s}{y 3 , y 4 , u} \ x l5 x 4 , t) denotes the probability density 
conditioned upon an event at later time t that the components for I = 1, 2 of the system 
(|H]) take the value yi,y2 at a prior time s whilst those for / = 3,4 the value y3,y 4 at 
time u < t. Since we are interpreting fllTf ) as BSDEs with final conditions this event is 
perfectly well defined. Therefore we conclude: 

d t «nf =1 T(x ; ,t) >/> v >= 

= -M 4 « Uf =1 T(x h t) > f > v > +C(x x -x 2 ) « T(x 3 ,t)T(x 4 ,t) > />v + 

+ perm. (28) 

In the limit t ] oo we obtain: 

M 4 « T(xi) T(x 2 ) T(x 3 ) T(x 4 ) > />v = 
= C(x! — x 2 ) << T(x 3 ) T(x 4 ) >/> v +permutations (29) 

The Hopf's identity satisfied by the general iV-points correlation function can be de- 
rived in the same way. 



5 Conclusion 

By means of the theory of stochastic differential equations the Hopf's identities have been 
proven for a passive scalar advected by a Gaussian delta-correlated in time velocity field 
for both the compressible and incompressible case. 
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